Brownian motors: current fluctuations and rectification efficiency 
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With this work we investigate an often neglected aspect of Brownian motor transport: The role of 
fluctuations of the noise-induced current and its consequences for the efficiency of rectifying noise. 
In doing so, we consider a Brownian inertial motor that is driven by an unbiased monochromatic, 
time-periodic force and thermal noise. Typically, we find that the asymptotic, time- and noise- 
averaged transport velocities are small, possessing rather broad velocity fluctuations. This implies 
a corresponding poor performance for the rectification power. However, for tailored profiles of the 
ratchet potential and appropriate drive parameters, we can identify a drastic enhancement of the 
rectification efficiency. This regime is marked by persistent, uni-directional motion of the Brownian 
motor with few back-turns, only. The corresponding asymmetric velocity distribution is then rather 
narrow, with a support that predominantly favors only one sign for the velocity. 

PACS numbers: 05.60.Cd, 05.40.-a, 05.45.-a 
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I. INTRODUCTION 

The channelling of particles by harvesting the ther- 
mal noise gives rise to a diffusive transport of particles. 
In periodic potentials the second law of thermodynam- 
ics implies that no net transport occurs. The situation 
changes drastically, however, in the presence of unbiased 
non-equilibrium noise which acts additionally on such 
systems. Then, the concept of Brownian motors does 
provide the possibility for directed, noise-induced trans- 
port. The phenomenon has widespread applications in 
physics, chemistry and in the biological sciences where it 
can be put to work for shuttling reliably and efficiently 
particles on the micro-scale, or even on the nano-scale 

The vast majority of works on Brownian motors con- 
centrate on the behavior and the selective control of the 
emerging directed transport as a function of control pa- 
rameters such as the temperature T, an external load F 
(yielding the load-current characteristics), or some other 
control variable. In contrast, the role of the fluctuations 
of the directed current has not attracted much attention 
in the literature. A notable exception is the first work on 
an inertial (rocking) ratchet ^ wherein the higher order, 
statistical cumulant properties of the stochastic position 
variable have been explored. Here, we fill this gap and 
focus in more detail on the fluctuating behavior of the 
Brownian motor current. The average drift motion to- 
gether with its fluctuation statistics are salient features 
when characterizing the modus operandi of a particular 
Brownian motor. 

It is intuitive that the fluctuations of the drift variable 
do impact the overall " efficiency" of the transport under 
consideration. The objective of an optimal operation of a 
Brownian motor machine can be formulated in a variety 
of ways, e.g. see in Q- In close analogy to heat-engine 
machines, one can define a generalized, nonequilibrium 
efficiency of a Brownian motor. There exists no univer- 



sally agreed upon definition of this notion 0, 0, S ~ 
for a recent review on efficiency of Brownian motors see 
Ref. 13 . The most common definition of the efficiency of 
a Brownian motor operating at non-equilibrium is based 
on the ratio of the work (or power) done by the parti- 
cle against an external load and the input energy (input 
power). With this working definition the load force is 
inevitably included; in particular, this yields the result 
that the efficiency assumes a zero value when no load 
force is acting. 

Alternative proposals for efficiencies have been pro- 
posed as well la, IE • Some of these proposals do pro- 
vide a nonzero value for a vanishing bias force. Yet 
another possibility consists in characterizing the recti- 
fication power of the transport in terms of the so-called 
Peclet number, i.e. the quotient of the drift velocity and 
the diffusion. This notion has been used in continuous 
state Brownian motor transport , and also for discrete 
motor models The concept is related in spirit to 

a Fano factor measure of the velocity fluctuations 
used, for example, to characterize molecular motors ^3 ■ 

In this work, we shall follow the reasoning of Suzuki 
and Munakata , in order to characterize the efficiency 
of rectification in absence of external bias forces. Typ- 
ically, there occurs a competition between two mecha- 
nisms: a "giant enhancement" of diffusion ITolIT^ [l5l | 
and an optimally large, (uni)-directional transport veloc- 
ity [l^ |l3| . The first perspective aims at controlling the 
magnitude of the effective diffusion independently of the 
temperature. It thus carries a rich potential for techno- 
logical separation devices. The second facet attempts to 
achieve a maximal "coherence" for the transport. Such a 
coherence is of relevance for Brownian motors modelling 
biophysical molecular motors ISj. 

Most of the Brownian motors and the majority of ap- 
plications studied in the prior literature operate in the so- 
called overdamped Brownian motion regime. For specific 
applications the role of inertial effects can become, how- 
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ever^f salient importance HIHITlliiilllillpil 
I25I l2a . I27I l2a | . The overdamped dynamics is a valid ap- 
proximation for many physical applications . It is also 
particularly well suited to describe the motion of molecu- 
lar motors In other situations, however, the iner- 
tial effects cannot be neglected. An exemplar is the diffu- 
sion of atoms on a crystal surface |29j . There, the dynam- 
ics may be underdamped, exhibiting long correlated hop- 
ping among binding sites. This physics has been verified 
experimentally by use of scanning tunnelling microscopy 
[sof , field ion microscopy [sj , or for quasi-clastic Helium 
atom scattering [s^. 

The inclusion of inertia adds significant complexity to 
the problem. This is so, because a periodically rocked, 
single degree of freedom with nonzero mass possesses a 
three-dimensional phase space that can exhibit a chaotic 
dynamics 0, |23, |^ . This is in contrast to the case of 
rocked, overdamped Brownian motors 's^, "s?, "sF, "sg", 'st'I . 
While the chaotic dynamics of a driven-damped particle 
in a symmetric periodic potential has been investigated 
thoroughly during the 1980's, the (chaos)-induced, di- 
rected transport of an asymmetric, inertial Brownian mo- 
tor has been pioneered only much later in There, it 
has been demonstrated that the corresponding dynam- 
ics features a rich structure, possessing many intrigu- 
ing current-reversals. In the deterministic case, different 
asymptotic solutions can coexist, e.g. running and locked 
states. Moreover, the onset of diffusive behavior due to 
chaotic dynamics has been investigated in terms of the 
second moment of the particle position diffusion. This 
first study was followed up with more detailed investi- 
gations ^u\, where it has been shown, for example, that 
the transport may reverse the direction at the transition 
from a chaotic to a regular motion. Additionally, inter- 
mittent trajectories have been observed. In such cases 
the system follows for a certain time a regular orbit, but 
then suddenly switches to a sticking orbit. The result- 
ing average flux depends on the time which the particle 
spends in a particular state, and the system may exhibit 
super-diffusive behavior. In related work '24'|, the inertial 
Brownian motor dynamics in a regime of adiabatic driv- 
ing has been investigated. These authors focused on the 
onset of the diffusive transport as the damping coefficient 
decreases. 

In this work we concentrate on the connection between 
the directed transport and its fluctuation characteris- 
tics. This study is of relevance for the optimization of 
rectification: The directed current should not become 
swamped with the unavoidable fluctuations of the trans- 
ported quantity. Our investigation is based on an inertial, 
noise-driven rocking ratchet. We shall analyze quanti- 
ties such as: the long-time averaged velocity, its velocity 
fluctuations, the fluctuations of the kinetic energy of the 
Brownian motor and the efficiency for rectification. 

The paper is organized as follows. In the next section 
we present the Brownian motor model. In Section III, we 
elaborate on the problem of the efficiency for rectifying 
noise in connection with the fluctuation behavior of the 



Brownian motor current. In Section IV, we describe our 
numerical flndings for a generic set of parameters, while 
in Section V, we elucidate the optimal working conditions 
for rectification and directed transport. 

II. MODEL 

To start, we consider the motion of a classical particle 
of mass m moving in the periodic, asymmetric ratchet 
potential V{x). The particle is driven by an unbiased 
time-periodic, monochromatic force of strength A and an 
angular frequency Vt. The dynamics is additionally sub- 
jected to thermal noise. The Brownian motor dynamics 
is thus governed by the Langevin equation [39l| 

mi + ^i^-V'{x)+Acos{nt) + ^/2^^{t), (1) 

where the prime denotes a differentiation with respect to 
the argument of V{x). The parameter 7 is the friction 
coefficient, T denotes temperature, and k is the Boltz- 
mami constant. The ratchet potential V{x) — V{x + L) 
has the period L and a barrier height AV. Thermal 
ffuctuations are modelled by the zero-mean (5-correlated 
Gaussian white noise £,{t). This noise term obeys the 
Einstein relation with the noise correlation given by 
{£.{t)£,{s)) — S{t — s). We next introduce dimension- 
less variables. The natural length scale is given by the 
period L of the ratchet potential. The dynamics pos- 
sesses several time scales. We introduce the character- 
istic time To as determined formally from the Newton 
equation, mx ~ —V'{x), by balancing the two forces 
mL/T§ = AV/L; yielding = niL^/AV. The scaled 
variables thus read: 



The dimensionless Langevin dynamics consequently as- 
sumes the form 

i + ji^ -V'{x) +acos{Lui) + ^/2jDo (3) 

where 

• the re-scaled friction coefficient 7 = {'~f/m)To is the 
ratio of the two characteristic time scales, tq and 
the relaxation time scale of the velocity degree of 
freedom, i.e., — m/7, 

• the potential V{x) — V{x)/AV assumes the period 
1, and the barrier height equals unity; i.e. AV^ = 1, 

• the drive has the re-scaled force strength a = 
AL/AV with the dimensionless angular frequency 

• the re-scaled, zero-mean Gaussian white noise 
forces obey {£,{i)^{s)) = 6{i — s) with a re- 
scaled noise intensity Dq = kT/AV. 
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In the following, mostly for the sake of simplicity, we shall 
only use dimensionless variables and shall omit the "hat"- 
notation in all quantities. For the asymmetric ratchet 
potential V{x) we consider a linear superposition of at 
least two (or even three) spatial harmonics (see in Sect. 
IV), i.e., 

V{x) = Vo[sin(27rx) + ci sin(47ra;) + C2 sin(67r2;)], (4) 

wherein Vq normalizes the barrier height to unity, and the 
parameters ci and C2 characterize the spatial asymmetry. 



III. FLUCTUATION AND RECTIFICATION 
MEASURES 

Throughout the following we focus on the asymp- 
totic, periodic regime after effects due to the influence 
of initial conditions and transient processes have quiet 
down. Then, the main statistical quantifiers of the driven 
stochastic process can be described in terms of time- 
and ensemble-averages. For a given quantity f{x{t)) 
its time-homogeneous statistical properties are obtained 
only in the long-time limit after transients have died out 
and after both, the average over the temporal period 
of the driving and the corresponding ensemble-average 
are performed [3^ . In this asymptotic regime the time- 
independent (single-time) quantities are obtained by a 
double averaging procedure over both the noise and the 
period of driving; i.e., 



lim — 



t+27r/w 



< f{x{s)) >- ds, 



(5) 



wherein -< /(s) >~ indicates the average over the noise 
realizations (ensemble-average) . 

The most salient transport quantity is the average, di- 
rected velocity (v) of the driven Brownian particle. Here 
v{t) denotes the stochastic process x{t) in Eq. 13 Of 
equal importance are, however, the fluctuations of v{t) 
around its mean (v) in the long time regime, i.e., the 
variance 



(6) 



The Brownian motor moves with current values v{t) that 
range typically within 



(7) 



If ay > (u), and even more so if a^, ^ (u), the Brownian 
motor can possibly move for some time in the opposite di- 
rection of its average value (v). The question thus arises: 
Is an efficient directed transport still feasible? 

To answer this challenge we shall introduce a measure 
for the efficiency 77 of the rectification of thermal noise, 
a quantity directly related to the velocity fluctuations. 
Here, we follow the reasoning of Suzuki and Munakata 



0, which yields a nonvanishing rectification efficiency 
also in the absence of an external bias. This efficiency 
of rectification follows from an energy balance of the un- 
derlying inertial Langevin dynamics. When specialized 
to our situation, rj is given by the ratio of the dissi- 
pated power 7(w)^ associated with the directed motion 
of the motor against friction, and the input power from 
the time-periodic forcing. The result assumes the explicit 
form (see in the Appendix) 



\{v^)-Do\ ■ 



(8) 



It thus follows that for a decreasing variance of the 
velocity fluctuations, a^, the efficiency of the Brownian 
motor increases. This is just what one would expect on 
naive grounds: The transport of a Brownian motor can 
be optimized in regimes of a large, directed average cur- 
rent which intrinsically does exhibit only small fluctua- 
tions. Moreover, in our study we found numerically that 
(w^) > Do holds true for any chosen set of the simulation 
parameters. 



IV. FLUCTUATION BEHAVIOR OF CURRENT 
IN AN INERTIAL ROCKED BROWNIAN 
MOTOR 

Deterministic inertial Brownian motors exhibit a com- 
plex dynamics including chaotic regimes. The applica- 
tion of noise then generally destroys the complex flne 
structure of their phase space and tends to smooth out 
their characteristic response function. 

There are two classes of states of the driven system 
dynamics: the locked states, in which the particle stays 
inside one potential well, and the running states for which 
the particle runs over the potential barriers. The first 
regime is characteristic for small driving strengths. When 
the amplitude of the external force is made sufficiently 
large we find that running states appear. These running 
states can either be chaotic (diffusive) or regular. 



A. Numerical schemes 

We have carried out extensive numerical studies in or- 
der to identify generic properties of the noise-induced 
transport. Applying a small, but finite noise strength, 
Dq > 0, we have integrated the Langevin equation (PJ 
by employing the Euler method with a time step of 
h = 10^^. For the initial condition of the coordinate 
x(t) we used a uniform distribution over the dimension- 
less period L = 1 of the ratchet potential. Likewise, 
the (scaled) starting velocity has been chosen at random 
from a symmetric, uniform distribution over the interval 
[—1, 1]. All quantities of interest were averaged over 250 
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different trajectories. Each single trajectory evolved over 
45 X 10"^ periods. The transient regime usually relaxed al- 
ready long before 500 periods of the driven dynamics had 
elapsed. For the cases of very weak noise and weak driv- 
ing we extended the corresponding time-span to ensure 
that the transient dynamics has quiet down completely. 
In the limiting deterministic case, i.e. Dq = 0, we used 
the Runge-Kutta algorithm of order 5. In this case, the 
averages were calculated over 10'^ differing trajectories, 
each trajectory evolving over 10'^ periods. 



B. Numerical results: Current fluctuations versus 
driving strength 

We start out to study the role of fluctuations by vary- 
ing the amplitude a of the sinusoidal driving force. In do- 
ing so, we assume a relatively small temperature, so that 
the Brownian motor dynamics is not far from a determin- 
istic behavior as described in prior works 2, 20, 24] . The 
average asymptotic long-time current velocity is shown 
in Fig. n^a). It reveals that for an amplitude a ~ 1.5 
the directed, inertial transport sets in before the lower 
threshold of the ratchet force is reached. It assumes a 
first local maximum near the lower threshold of the po- 
tential force near a ~ 2.1. Note that in the presence of 
small noise, the current is strictly speaking never zero. 
For all practical purposes, however, we can characterize 
the outcome of our Langevin simulation as a determinis- 
tic, zero-current result. Below this threshold the system 
mainly dwells in the locked state. Upon closer inspection, 
we notice that in the vicinity of a ~ 0.6, the velocity fluc- 
tuations (Ji, shown in Fig. ^b) undergo a rapid increase. 
In Fig. nic) we also display the relative fluctuations of 
the kinetic energy, E = v'^ /2, (the re-scaled mass is one); 
i.e. 

Around this value of the driving amplitude, such a quan- 
tity undergoes a giant enhancement. Finally, we remark 
that for the equilibrium Maxwell distribution we flnd 
a^KE)"^ = 2, as expected. 

Upon further increasing the amplitude of driving, a > 
1.5, the Brownian motor generates for this set of pa- 
rameters the desired, directed transport behavior. At 
the same time we observe that the width of the weakly 
asymmetric, corresponding distribution P{v) slightly de- 
creases, meaning that the velocity fluctuations become 
smaller. The following explanation thus applies: Because 
at a < 1.5 escape jumps between the neighboring wells 
are rare, i.e., the average directed current is very small 
(note also the accompanying, very weak asymmetry in 
the velocity distribution). The input energy is pumped 
primarily into the kinetic energy of the intra-well mo- 
tion and eventually dissipated. As a is increased further, 
the Brownian motor mechanism starts to work, and some 
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FIG. 1: (Color online) Fluctuation behavior of an inertial 
Brownian motor versus the driving strength a. (a): av- 
eraged dimensionless velocity {v) of the inertial Brownian 
motor in Eq. (|3J; (b); variance of the corresponding ve- 
locity fluctuations (t„; (c): fluctuations of the rescaled ki- 
netic energy a%/{E)^\ (d): rectification efficiency in Eq. ||5J. 
All quantities have been computed for the rescaled poten- 
tial V{x) = -Vo[sin(27ra;) + 0.25 sin(47r2:)], where Vo ~ 0.454 
normalizes the barrier height to unity, see the inset in (a). 
The force corresponding to this potential ranges from —2.14 
to 4.28. The angular frequencies at the well-bottom and at 
the barrier-top, respectively, equal each other, reading 5.28. 
Bottom panel (e): velocity distributions for selected driving 
amplitudes, i.e. a = 0.68, 1.5, 2.14, 3.25. All these distribu- 
tions were normalized by setting their maximum to one. The 
remaining rescaled parameters read: friction 7 = 0.5, angular 
driving frequency u = and weak thermal noise of strength 
Do = 0.01. 



part of energy contributes to the net motion of the par- 
ticle. Therefore, less energy remains available to drive 
intra-well oscillations and consequently the distribution 
P{v) shrinks, see in Fig. ^e) 

Correspondingly, due to inertia, the mean velocity in- 
creases, reaching a second maximum before the upper 
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threshold value of the potential force a ~ 4.28. Above 
this driving amplitude, the current starts to decrease be- 
cause of the weakening influence of the ratchet potential 
at large rocking amplitudes. 

The occurrence of multiple reversals for the directed 
current, as it occurs in Fig. ^a), is a known, interesting 
feature of inertial Brownian motors. Several prior studies 
did elucidate in greater detail the corresponding mecha- 
nism at work |2. 120ll22. 1241 l25l | . Here, we take instead a 
closer look at the current fluctuations. We observe that 
for the chosen set of parameters the maximal stationary 
velocity in Fig. QJa) does not exceed the value 0.4. In 
contrast, its fluctuations keep growing as the driving am- 
plitude rises. At large driving, the particle no longer feels 
the potential and undergoes a rocked, free Brownian mo- 
tion with the velocity fluctuations growing proportional 
to a, cf. Fig. mb). 

On the other hand, the relative fluctuations of the ki- 
netic energy do saturate, see Fig. ^c). These are sup- 
pressed to values near 1, which lies below the equilibrium 
value of 2. Within this directed transport regime, the ef- 
ficiency © remains rather small, cf Fig. QJd). Such small 
rectification efficiency is the rule for this driven inertial 
Brownian motor. 

Let us next inspect the current probability distribution 
P{v). All P{v) curves reported in the following have 
been normalized so that their maximum (i.e., their high- 
est peak) is set to a fixed, unit value. Only then we 
can detect the details in their shape upon varying the 
corresponding parameter such as the driving amplitude 
a or the noise strength Dq- These probabilities look 
rather symmetric; however, a finite ratchet velocity re- 
quires a certain amount of asymmetry either in the lo- 
cation and/or the width of the velocity peaks. Here, the 
current results mainly due to a slight shift of the maxima 
location. 

The most peculiar feature of the current distributions 
shown in Fig. ^^e) is the emergence for a > 0.6 of two 
additional side-peaks centered near w = ±1, which even- 
tually dominate P{v) at larger driving amplitudes. Of 
course, for zero drive P{v) boils down to a single-peaked 
Maxwell distribution, strictly symmetric around v = 0. 
To investigate the onset of these two side-peaks for a val- 
ues corresponding to vanishingly small currents, we set 
P{v) = qPo{v - 1) + qPo{v + 1) + (1 - 2q)Poiv), where 
q varies from to 1. If Pq is taken to be a symmetric 
Gaussian function, then the kinetic energy fluctuations 
can be evaluated; ctb exhibits the behavior shown in Fig. 

CTc). 

What it the origin of those three peaks in the distri- 
bution P{v)7 Our first conjecture to connect it with the 
'running' solutions turned out to be incorrect. This is so, 
because for a 1 the particle rarely leaves the confining 
potential well and thus cannot significantly contribute 
to the side peaks of the distribution function. We fur- 
ther checked the outcome for the velocity distribution 
when reflecting barriers were placed at the maxima of 
the potential. Under such constraints, the three-peak- 



structure is recovered as well. Moreover, the sinusoidally 
driven damped particle in a harmonic potential can ex- 
hibit both, a singly-peaked as well a doubly-peaked av- 
eraged velocity distribution, see in Ref. [38j. However, 
for the parabolic potential that fits best the wells of our 
ratchet potential around its minima, we found a single 
peaked P{v). 

We therefore do conclude that the characteristic be- 
havior for the additional side-peaks is rooted in the non- 
linear, anharmonic character of the corresponding well of 
the periodic asymmetric ratchet profile. 



C. Numerical results: role of the shape of the 
underlying ratchet potential 

We next study the influence of the ratchet profile (or 
force) for the mean velocity and the corresponding fluctu- 
ation behavior of the directed current. We use a stylized 
potential shape composed of three spatial higher harmon- 
ics with C2 ^ 0, see in Q. Note that this potential 
shape possesses an opposite polarity as compared with 
the ratchet potential depicted in the inset of Fig. ^a). 
Put differently, the natural direction of the Brownian mo- 
tor motion is to the left, in the direction of the weaker 
slope, see Fig. |2Ia), inset. 

For small driving amplitudes, this inertial motor pre- 
dominantly dwells in a potential well. The directed cur- 
rent is very small and negative (directed towards the left). 
If the driving amplitude exceeds the upper threshold am- 
plitude of the ratchet force at a ~ 4.66, the motor starts 
to move more regularly, reaching an extremal speed near 
a ~ 6. Interestingly enough, the motor moves now to- 
wards positive a;-values. The resulting transport velocity 
thus cannot be easily predicted a priori in this nonadi- 
abatic driving regime. This is a benchmark feature of 
these inertial Brownian motors where the coupling be- 
tween the deterministic driving force, F(t) = acos{Lut), 
and the resulting motion of the driven Brownian particle 
are coupled loosely, only. Except for a narrow regime of 
sizable values, cf. Fig. |2ta), the emerging average ve- 
locities are typically very small, yielding no correspond- 
ing efficient rectification. Therefore, the efficiency of this 
more complex Brownian motor mimics again closely the 
behavior of the average motor velocity, see panel (d) in 
Fig.ia 

The velocity fluctuations exhibit a similar behavior as 
in the case of Fig. ^discussed above. The variance grows 
nearly linearly with increasing driving strength. Typi- 
cally, the average velocities are small and the fluctuation 
behavior is similar to the behavior discussed above for 
the first ratchet potential. Interestingly enough, however, 
above threshold driving induces a distinct peak behavior 
for the velocity which is accompanied by a corresponding 
dip in the behavior for the velocity fluctuations. This dip 
in the variance then gives rise to a window with appre- 
ciable efficiency, cf. Fig. Hfd). 

The behavior for the velocity distributions is again 
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FIG. 2: (Color online) Tailoring the shape of the potential, 
(a): averaged dimensionless velocity {v) of the inertial Brow- 
nian motor under nonadiabatic driving conditions; (b): cor- 
responding velocity fluctuations a^; (c): fluctuations of the 
rescaled kinetic energy a%/{E)'^; (d): efficiency. All quan- 
tities are plotted versus the external driving amplitude a 
for the asymmetric ratchet potential V{x) = Vo[sin(27ra::) + 
0.245 sin(47rx) -I- 0.04 sin(67ra;)], where Vb — 0.46f normalizes 
the barrier height to unity (see inset in (a)). The forces stem- 
ming from such a potential range between —4.67 and f.83. 
The two angular frequencies at the well-bottom and at the 
barrier-top are the same, reading 5.34. The corresponding 
velocity distributions P{v) axe displayed in panel (e) for the 
indicated driving amplitudes, i.e., a = 0, 4, 6, 7. All P{v) 
curves have been normalized by setting their maximum to 
one. The remaining parameters are: 7 = 0.9, uj = 4.9 and 
Do = 0.01. 



generic: Small average velocities exhibit nearly symmet- 
ric velocity distributions, see Fig. Ole). Only for the 
nonadiabatic peak behavior of the mean velocity does 
one identify also an appreciable asymmetry for the ve- 
locity distribution. 



D. Numerical results: Current fluctuations versus 
noise strength Dq 

In Fig. 13 we numerically investigate the directed 
transport versus the temperature Dq. We have chosen 
a sub-threshold driving stren gth for which the thermal 
noise plays a constructive role [4(j| by inducing noise acti- 
vated jumps across the potential barriers. We set a = 0.8 
and the other parameters remain the same as in Fig. ^ 
Then, we find a characteristic velocity reversal near the 
dimensionless temperature Dq = kT / IS.V ~ 1. There, 
the thermal energy compares with the activation energy 
over the barrier height of the ratchet potential. A subse- 
quent increase of temperature causes a diminishing role 
of the asymmetric ratchet potential and, consequently, 
the directed transport degrades. 

Moreover, the time-averaged velocity distribution ap- 
proaches the equilibrium velocity distribution Q. It is 
remarkable that within a certain range of temperatures 
the fluctuations of the kinetic energy exceed the rele- 
vant value for the Maxwellian equilibrium distribution, 
cf. Fig. ini^c). A shallow, local minimum occurs for 
the velocity fluctuations where the average current it- 
self is maximal. These fluctuations are, however, notably 
three orders of magnitude larger than the small-valued, 
directed current. Not surprisingly, the rectification ef- 
ficiency shown in Fig. |3fd) is quite small. Again, the 
Brownian motor is not operating optimally. 



V. TAILORING RECTIFICATION EFFICIENCY 

Thus far, changing the ratchet profiles did not lead to 
a large enhancement of the rectification efficiency. What 
is needed in achieving a large rectification efficiency is 
a sizable Brownian motor current which is accompanied 
by small current fluctuations only, see Eq. ((SJ. This sce- 
nario seemingly implies that the directed current should 
proceed in a persistent manner with very few, occasional 
back-turns only. This in turn causes small fluctuations 
in the velocity and, additionally, provides a dominating 
asymmetry of the velocity distribution. 

Such a behavior can be realized by a combined tailor- 
ing of the asymmetry of the ratchet potential together 
with the use of appropriate driving conditions. In the 
quest for achieving such a favorable situation we use the 
three-harmonics ratchet potential plotted in the inset of 
Fig. |2Ia). Our hope is that upon minimizing the noise 
further we can achieve a substantial improvement of the 
efficiency. 

At very weak noise and large, nonadiabatic rocking 
frequencies, this inertial Brownian motor starts moving 
efficiently near the upper threshold of the ratchet force 
a ~ 4.66, see Fig0] (a). Because the directed velocity 
becomes maximal and simultaneously its fluctuations are 
locally minimal, see in Fig. EJb), we indeed flnd the 
desired enhancement of the rectification efficiency, see 
Fig. Efd). The fluctuations of the kinetic energy grow 
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FIG. 3: (Color online) Fluctuation behavior of an inertial 
Brownian motor versus the noise strength Dq. (a): averaged 
dimensionless velocity (i;); (b): the corresponding velocity 
variance CTi,; (c): fluctuations of the rescaled kinetic energy 
a%/{E)'^; (d): the corresponding efficiency. Numerical results 
obtained for the same ratchet potential as in Fig. [H In the 
bottom panel (e), the velocity distribution P{v) is shown for 
Do = 0.01, 0.1, 2, 10. All P{v) curves are normalized as in 
Fig. The remaining rescaled parameters are: 7 — 0.5, 
uj = 3.6 and a = 0.8. 



slightly; nevertheless, these are still strongly suppressed 
in comparison to the equilibrium value 2. 

We have studied several other ratchet potentials by 
varying the parameters ci and C2 in Eq. Q and still found 
regimes where the inertial ratchet works with a high effi- 
ciency (not shown). In all these cases we found that the 
velocity distribution has a support concentrated mainly 
on one of the semi-axes. Strongly asymmetric velocity 
distributions are depicted with Fig. EJe). In contrast, 
with the mode C2 set zero (see in inset in ^ we could 



not identify such an optimal regime for rectification of 
noise. The shape of these distributions just corroborates 
the fact that large rectification efficiencies are the result 
of persistent, (uni)-directional Brownian motor motion, 
accompanied by a strong asymmetry of the current statis- 
tics. 



VI. CONCLUSION 

With this work we have elucidated, directed Brownian 
motor transport in rocked rachet potentials in the pres- 
ence of inertia and thermal noise. We focussed on several 
parameter regimes and studied by numerical means the 
operation of this massive ratchet machine. In particu- 
lar, we investigated the variation of the average current 
versus driving amplitude a and the temperature strength 
Dq. Our main objective has been the behavior of the ac- 
companying current fluctuations as a function of these 
transport parameters. These fluctuations crucially im- 
pact the rectification behavior, as measured by the rec- 
tification efficiency in Eq. ((SJ. 

Typically, the current values and the corresponding ve- 
locity fluctuations are so, that no appreciable rectifica- 
tion emerges in these inertial, rocked Brownian motors. 
There exist, however, tailored regimes of rachet profiles 
and driving parameters for which an enhancement of rec- 
tification and optimal transport do occur. These regimes 
are marked by a large Brownian motion transport with 
few back-turns only. This in turn implies a narrow, asym- 
metric velocity statistics with dominantly, one-sided sup- 
port of either positive- or negative- valued velocities. 

These novel findings for the fluctuation statistics of 
Brownian motor velocities can be put to use in diverse 
technological devices that pump and separate efficiently 
and reliably Brownian particles in correspondent physical 
and biological Brownian motor systems 18]. More- 
over, the results derived herein for such driven inertial 
Brownian motors can be applied as well to the phe- 
nomenon of Stochastic Resonance [40l | in corresponding 
underdamped regimes. 

Appendix 

With this appendix we present the derivation of the 
expression ((SJ for the rectification efficiency. We combine 
the arguments in [sj, ^ 0| and establish the efficiency ?/ 
as 



In the denominator. Pin denotes the rate of the energy 
input to the system. There is no overall consensus on the 
numerator ^ 0, 0, IE • HA denotes the rate of work 
done on the fluid by the Brownian motor motion, then the 
corresponding efficiency rj is not an appropriate measure 
because A {v^)- This quantity can be relatively large 
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FIG. 4: (Color online) Tailoring the shape of the potential, 
(a): averaged dimensionless velocity {v) of the inertial Brow- 
nian motor under nonadiabatic driving conditions; (b): cor- 
responding velocity fluctuations av; (c): fluctuations of the 
rescaled kinetic energy a%/{E)'^; (d): efficiency. All quanti- 
ties are plotted versus the driving amplitude a for the asym- 
metric ratchet potential of Fig. |2| The corresponding velocity 
distributions P{v) are shown in panel (e) for selected driving 
amplitudes, i.e., a = 0, 3.25, 3.5, 4.5. All P{v) curves are nor- 
malized as in Fig. |5| The remaining parameters are: 7 — 0.9, 
Lu = 4.9 and Do = 0.001. 



where V{x,t) = V{x)—ax cos{ujt) andW{t) is the Wiener 
process {{W{t)} = 0, {W^{t)} = t). 

Now, we evaluate the ensemble and temporal averages 
of the re-scaled kinetic energy G{v) = v = v{t). To 

this aim, first we apply Itos differential calculus to the 
function G{v) to obtain 

d[v'^/2^ ^ -(^jv^ +vV'ix,t) -"fDa^dt (13) 
+ ^/2-fDovdW{t). 

The ensemble average (i.e. the average over all realization 
of the Wiener process denoted by -< ... y) for the rate of 
change of the kinetic energy results in 



dt 



- ^ ^ v'^ y + ^ vV'{x) y- 
< V acos{Lut) > jDq 



(14) 



where we exploited the (Ito)-martingale property (for 
the part containing the Wiener process). Next, we av- 
erage over the temporal period as in (O (periodic time- 
dependence of asymptotic probability). In doing so, we 
evaluate 



{—v^)^<v'^{t + 2T:/Lu)> <v^{t)'^=0. (15) 



Likewise, for the contribution 

{vV'{x)) =< V {x{t + 271 /u)) > <V {x{t)) 

Consequently, we obtain 



(16) 
(17) 



where the combined average Pin = {v{t)aco8ujt)) is the 
input energy to the system per unit time. Thus, upon 
combining (|10|l and (|17() the relation in (jSJ emerges. We 
also emphasize here, that our scheme for the efficiency of 
rectification at zero bias is independent of the transport 
friction-coefficient 7. This feature is in agreement with 
the corresponding result by Suzuki and Munakata j3j . 
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